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We give a detailed study of axion-photon and photon-axion conversion amplitudes, which enter 
QQ . the analysis of "light shining through a wall" experiments. Several different calculational methods 

' are employed and compared, and in all cases we retain a nonzero axion mass. To leading order, we 

, find that when the photon frequency u) is very close to the axion mass m, there is a threshold cusp 

^SJ ■ which significantly enhances the photon to axion conversion amplitude, by a factor tj/Vi^^ — m-^ 

relative to the corresponding axion to photon conversion process. When m = Q, the enhancement 
I factor reduces to unity and the results of previous calculations are recovered. Our calculations 

include an exact wave matching analysis, which shows how unitarity is maintained near threshold 
at Lj = m, and a discussion of the case when the magnetic field extends into the "wall" region. 

in 

Dh; I. INTRODUCTION 

D ' The axion is a remarkable idea proposed at the end of the seventies in order to solve the strong CP problem 
1-^ ' The possible existence of the axion could help solve long-standing cosmological problems and, therefore, searching for 
it is an important issue 13]. 

. Following on the seminal papers of Sikivie Q which showed that axions can be detected through axion-photon 
J> ■ conversion in a magnetic field, twenty years ago Van Bibber et al. proposed a "light shining through a wall" (LSW) 
, experiment where photon regeneration could be used as an indication of an axion-photon coupling [3| . In the Van 
' Bibber et al. setup (see Figure 1), B is an external magnetic field and the wall is opaque for photons while it is 
transparent for the very weakly interacting axions. An x-ray version of this experiment has recently been proposed 
by Rabadan et al. 5] , and two recent optical laser LSW experiments have reported preliminary results 1 H • 

Motivated by the current interest in LSW experiments and axion detection, we give here a detailed examination 
of the theory of photon-axion and axion-photon conversion processes. Our aims are in part to compare different 
calculational methods, and in part to examine threshold effects that appear when the axion mass is not set equal to 
zero. 

In the LSW experiments one can distinguish seven zones, as sketched the figure. We are mainly interested in the 
transition zones I-II-III and V- VI- VII, namely, vacuum-magnetic field-vacuum regions where an incident photon gives 
5J] ■ rise to an exiting axion, and an incident axion gives rise to a regenerated photon. Although we will mainly focus on 
\ the theory connected with these regions, we will also briefly discuss what happens when the magnetic field enters the 
absorbing wall region. 

The paper is organized as follows. In Sec. II we give the Lagrangian, equations of motion, and some elementary 
properties following from the equations of motion, such as the continuity conditions, free space kinematics, conserved 
probability current, and unitarity relation. In Sec. Ill we derive the leading order photon-axion and axion-photon 
conversion amplitudes by a Green function method, and in Sec. IV we repeat this derivation by a WKB/eikonal 
method, in both cases allowing the magnetic field to have a general dependence on z, the coordinate along the axis 
of the experiment. We find that when the axion mass is taken into account, the ratio of the photon-axion conversion 
amplitude to the axion-photon conversion amplitude is w/fc, with w and k = \fuP^—rr? respectively the photon and 
axion wave numbers. Since the leading order axion-photon amplitude violates unitarity near threshold at a; = to, in 
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Sec. V we do an all orders calculation for the case of a piecewise constant magnetic field. This gives the shape of the 
photon-axion conversion amplitude near threshold, and allows us to check that the unitarity constraint is obeyed. In 
Sec. VI we give further results following from the calculation of Sec. V. We show that the ratio of photon-axion to 
axion-photon amplitudes derived to lowest order is in fact exact at all orders for a piecewise constant magnetic field, 
we show that the all orders calculation restores unitarity near threshold, and we give some rough estimates of "light 
through a wall" and photon-axion conversion probabilities integrated over the threshold cusp region. In Sec. VII we 
briefly discuss what happens when the magnetic field penetrates the wall region; although the wall strongly absorbs 
photons, the weak coupling of photons to axions leads to a greatly suppressed absorption of the axion wave. Finally, 
in Sec. VIII we give a brief summary and suggestions for follow-up work. 



II. LAGRANGIAN, EQUATIONS OF MOTION, AND UNITARITY RELATIONS 

In this paper we will consider a system described by the action 



5 d'^x 



(1) 



with F'^" = Faf} , where </> is a real scalar field representing the axion, F^i, is the electromagnetic field tensor, 

and g is the coupling constant. 

The physical situation that we are interested in consists of propagating photons and axions in the presence of a 
static magnetic field as a background. We will suppose that this magnetic background B is pointing in the x direction, 
so that the relevant part of the interaction term is 

' d^xl3(l)E,, , (2) 

where /3 = gB and Ex is the component of the photon electric field parallel to the background magnetic field. 

Neglecting components of the photon field that do not couple to the axion through the magnetic field, and using 
Coulomb gauge, ([T]) becomes 

S = -\j [ad^a + (j){d'^ + m^)(j} - 2(3(j}dta] , (3) 

where a is a real potential which determines the relevant electric field component through E^ — dta. From this action 
we get the equations of motion, 

{d'^ + m'^)(P = Pdta , 

d^a = -~pdt<t) , (4) 

so that, specializing to the case in which the propagation direction is the z axis, we have 

{d^t ^dl + m^)4){z,t) = p{z)dta{z,t) , 

[d^^ - dl)a{z,t) = -P{z)dt^{z,t) . (5) 

The entire analysis of this paper will follow from the equations of motion ([S]). In particular, since these are second 
order in z, the functions 4'{z^t) and a(z,i), and their first spatial derivatives, must be continuous. Since no complex 
numbers enter the equations, we can follow the usual method of generalizing the wave amplitudes from real to complex. 
Multiplying the equation for by 0*, multiplying the equation for (jf by <j) and subtracting; similarly multiplying 
the equation for a by a*, multiplying the equation for a* by a and subtracting; and finally adding the two resultant 
equations, one finds a current conservation equation 

dtjo + dzjz = , 

jo — (p*dt(l) — dtcj)* (f) + a* dta — dta* a + I3{z){a* (j) — (f)* a) , 

jz = dz(j)*(t) — 4>*dz(j) + dza*a — a*dza . (6) 

Specializing now to the case of an incident wave with time dependence /(z)e^*'^*, where / denotes the fields <j> or A, 
jo becomes independent of time, and so dtjo — 0. The current conservation equation then becomes dzjz = 0, which 
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we use to give a useful unitarity relation, as follows. In free space regions where the magnetic field vanishes (zones I, 
III, V, and VII of Figure 1), the z dependence of A is e**'^^ and of <^ is e^*'"'^, where 

fc \/ up- — . (7) 

Assuming no incident waves from the right of the region where the B field is nonvanishing, the fields to the left will 
be given by 

a(z, t) = e"*'^*(a/e''^" + ai?e~''^") , (8) 
and the fields to the right of the B field region will be given by 

a(z, i) = e-*"*aTe''^^ . (9) 

Substituting ([8|) and ^ into ^ to give jz on the left and right, respectively, and equating these, we get the unitarity 
relation between incident, reflected, and transmitted wave amplitudes, 

+ uj\ai\^ = + \,^T?) + c^doflp + IflTp) . (10) 

This equation will play a role in our subsequent discussion. 



III. GREEN FUNCTION CALCULATION OF PHOTON-AXION AND AXION-PHOTON 

CONVERSION 

We turn now to a calculation of the photon-axion and axion-photon conversion amplitudes, to first order in the 
magnetic field strength parameter /3(z). 

By performing a perturbation expansion in powers of /3 we write the fields as follows 

0(z,i) = 0("'(z,i) + 0(i)(z,i) + --- , (11a) 

o(z,t) = a(°)(z,i) + a(i)(z,i) + --- , (lib) 

where 

(at'-<92 + m2)0(°)(z,t) = , (12a) 

(a2-c>2)a(o)(^^t) = , (12b) 

{dl~dl+'n?)<\}-^\z,t)- (i{z)dta!-'^\z,t) = , (12c) 
(52-a2)a(i)(z,t)+/3(z)at0(°)(z,t) = . (12d) 

In the last two equations, unperturbed fields [solutions of (|12ap and (jl2bp ] are sources for perturbed ones, so we 
will calculate the Green function for the differential operator in (|12a|l . Clearly, for case of the differential operator in 
(|12b[l it is enough to put to = 0. Thus, we are interested in the solutions Gm{z, t) of 

(d^-dl + m^)G^{z,t)^S{z)d{t) . (13) 

A direct calculation in Fourier space shows that the desired Green function is 

I r ^i(kz-Ljt) 

G^{z,t) = -j—^ dkdu- ^ , (14) 

[2itY J (w-tjfe)(w + Wfe) 



with ojk = -s/fc^ + m? ■ 

Since we want the retarded solution, that is Gm{z, t) = Q for t < 0, we integrate on the real w axis, with poles =F<^fe 
shifted to q^tj/c — ie. When i < 0, closing the contour up gives 0; when t > 0, closing the contour down gives the 
retarded Green function 



G'^{z,t) = -e{t) [—] I - U'ikz+'-.t) ^ e'ik^-^.t)\ ^ _oit)A^iz,t) . (15) 
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Writing A„i in a more compact form we obtain 



2-K UJk 



(16) 



and for the photon, taking to = 0, we find 



G'Jiz,t)^~9{t)Ao{z,t)^-0{t~\z\ 



(17) 



Returning now to the sohition of p2cp . for the case 0^°^ = of no incident axion, we take the incident photon field 
as a monochromatic wave a(°) = e^'^^^"*), so that the inhomogeneous equation (|12cP reads 



(df - + m^) 0(1) (z, i) = -itoPiz) e*^("-*) 
The sohition to p8|) . constructed using the Green function (|15p is 

rL 



(18) 



0(1) (z,i) = -tu! I dz'P{z')e'^'' 

rL 



l\—iu}t' 



iQe-"^' I dz'(3{z')e"^' 



dt'GZiz- z' ,t~t')e 



-oo Wfe V47r / 



(19) 



where we have assumed that (3{z) = outside the interval < z < L. 

The integral on t' can be done straightforwardly by defining T ~ t — t' and noting that, due to the 9 function, there 
are contributions only for T > 0. Then 



>{z,t) = -iwe-*"* / dz'(3{z')e"^' 



dz'p{z')e'^'' 



°° dk 

-oo ^fc V 47r 



2iu!k 



^ik{z — z') 



(20) 



The integral on k can be performed as usual in the complex plane. The boundary condition that a(°) vanishes in 
the infinite past can be included by giving uj a small positive imaginary part w ^ w + ie, e > 0, so that 



± \/ [uj + ieY ~ 'Ti^ ~ ± \/ llP- — rn? ± — = ± \/ uP- — vn? ± le' , e' > 



(21) 



We see that for z — z' > the contour can be closed up, circling the pole at A: = _ ra^ jg'^ while for z — z' < 
the contour can be closed down, circling the pole at A; = — if' . Equation ()20p thus gives 



)(z,t) 



2\/llP- — vn?- 



dz'e{z - z')f3{z')e 



UJ 



2Vi 



W — TO 





^-^Ot+^z^a--m- j dz' I3{z')e''' ^^^^^'^ 
L 







dz'/3(z')e*^'(^+^"'-'"') 



We see that the axion amplitude has a transmitted and a reflected part, 

0W(0,i) = 0SL(^,O + '/'RiLct(;2,t) , 

with 



(22) 



(23) 



2v;^ 



0Rilect(2;,t) 



2\/ llP- — 



icjt-\-iz\^uj^ — m'^ I p^^/^^iz' {uj—\^uj^ — m,^) 

Jo 

rL 



(24) 
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for general z, or for z values outside the magnetic field region, 



Lj2 — rn^ 

'/'iiLt(^,t) ^ ^ , ^ , / dz'(i{z')e'^'<'-^+^^^^^\ z<0 . (25) 

zyti^ — in Jo 

From psp we see that if the exponential in the integral for i^TrLa is rapidly varying on the scale over which f3{z) varies, 
the axion transmission is strongly suppressed. This statement can be given a quantitative form by an integration by 
parts, 

^ dz'Piz')e"'^''-^^^^'> r dz' ^ '}^'^ ^zz'(u-v^-~m-) ^ (26) 

where we have used the fact that /3(0) = /3(i) = to drop the surface terms. If we assume that \(3' [z) / (3{z)\ < 1/h, 
with h a measure of the characteristic distance over which the magnetic field varies, and that |/3(z)| < /3m, with /3m 
a measure of the maximum magnetic field, then substitution of (|26p into (|25p gives the inequality 

I^^Ll < ^ . " . (27) 

2 \l iJ^ — m?- h(ix) — \j llP- — rnP) 

Let us now consider /3(z) = /3, a constant. The previous integrals can be done straightforwardly, and we obtain for 
the transmitted part 



ct>i'L{z,t) = " gV(^-V5CT)(- ^ ^-2 _ ^2 ) -^Qj^^m^)] , (28) 

which has the modulus squared 



|41,(^, OP = 4 f -2 2 ) ('^ + - )' ^i^' - - )) ■ (29) 

Let us recall now that the unitarity relation of pop requires, for an incident photon amplitude of magnitude unity, 
that the transmitted axion amplitude obey the unitarity constraint 1 > (fc/a;)|(/)Tp, which in the notation of this 
section is 1 > {^/ui'^ — m? ■/w)|0(i)(z > L,i)|2. From ^ we see, however, that near threshold at w ~ r/i, we have 

{^uj-^-m?/uj)\(t)^^\z>L,t)\^ = sin'(mL/2) , (30) 



which approaches infinity ds, Q ^ m. Thus the lowest order expression for the transmitted axion amplitude violates 
unitarity near threshold; an all orders calculation, given below, is needed to see how unitarity is restored, and to give 
the shape of the axion transmission probability near threshold where k — \JuP' — rn? ~ 0^ /m. 

To complete the calculation, we will give the corresponding result for the photon field amplitude induced by an 
axion wave incident on the magnetic field region, with no incident photon (i.e., a^'^^ = 0). Using the zero mass 

retarded Green function of ^ with the source term -/3(z)<9f0(°)(z, t) = w/3(z)e*(^'^'"'"'^"'^*) in (fm]) . we have 







a(i)(z,t) = / dz' dt'G^iz - z',t- i')*'-'/3(-2')e'^^'^'""'^'"'^*'^ 



dz' / dt'-e{t -t' -\z- z'\) i<I./3(z')e*(^'^'-"'^'-'^*') 



J -OQ 

L poo 



dz' (3{z')e'' ""^^^ e^'^'^dT . (31) 

J\z-z'\ 

Finally, performing the integral onT — t ~ t' , we find for the transmitted photon amplitude dX z > 

ai]L{z,t) = -—— / dz'/3(z')e*"'(^"'-'"'-'=) . (32) 

2 Jo 

We see that when the axion mass is nonzero, the ratio of transmitted wave amplitude magnitudes is |(/'TrLs/flTrLs| — 
w/Vw^ — > 1. When the axion mass m vanishes, the transmitted amplitudes (|28p and ([5^ are identical up to a 
sign, as analyzed in greater detail by Guendelman [9] using his demonstration of an of an axion-photon duality. 
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IV. WKB/EIKONAL CALCULATION OF PHOTON- AXION AND AXION-PHOTON CONVERSION 

In this section we calculate the photon-axion and axion-photon conversion amplitudes, using a WKB/eikonal 
method. As an alternative method of solving the equations of motion, with e~"^* time dependence, we consider the 
following Ansatz 

cj}{z, t) = (j)o e'^(^)-*"*, a{z, t) = ao e^^^^)"*"* . (33) 

The equations of motion ([5]) now read 

{uj^ -m^ - + i0")<l>-ip(z)uja = , 

{oj^ ~ {x'f + tx")a + tP{z)uj<P = . (34) 

Expanding 9{z) and xiz) in powers of /3, 



e{z) = ^/uj^-m^z + e^^\z) + -- - , xiz)^ujz + x^^\z) + -- - , (35) 

where the index (1) indicates the first order term in the expansion in /?, equations (|34p become 



, (36) 



where for the moment we have dropped second derivatives of phases (as indicated below by a subscript WKB), and 
only first order contributions in p are considered. 

Non-trivial solutions satisfy the null determinant condition 



(eW)'(x«)' = T^C^ 

From the second row of (|36|) we get 



2 2 - (37) 



Ay/uj'^ — m 



= dz'(3{z')e'''^^^^^-^'> , (38) 

2ao J 



while from the first row of (|36p (or using the null determinant condition) we get 



0ii)^,)^Zl^ ^ / rf^.^(^.)g..'(--V^^) . (39) 

2<Po \J up- ~ vn? J 



This integral is the same as the one appearing in the transmitted axion wave in (I24p . In fact, we can now explicitly 
write the axion field to first order in /3, in the approximation of neglecting second derivatives in p6p . as 



V 200 - m'^ J J 



(40) 



or an equivalent expression, 



0„KB(z,t)-0(°n^,i)^ , , , g.(.y^^--t) / d//3(/)e-'(-^^^) . (41) 

The photon field satisfies a similar relation, namely 

a^^^{z,t)-a(°'>{z,t)~ _^e'-(-*) J' dz'P{z')e"'^^^^^^-^') . (42) 

With the substitution w — > cD, and taking the lower limit of integration as z = for the case when (3{z) vanishes for 
z < 0, these agree with the transmission amplitudes obtained by the Green function method. 



7 



Finally, let us perform the calculation, again up to first order in /?, but preserving the second derivatives of phases, so 
that the exact first order amplitude is obtained, including the reflection amplitudes calculated by the Green function 
method. From (l34l) we obtain 



which implies 



with the integral 



iMd'-^^y = iCe-^''^^"-'"'" +iaou;e- 



-lizyu. 



dz'f3{z')e 



(43) 



(44) 



(45) 



Integrating again, the phase O'^^^ is given by 

iC 



i^oO^^^ = C + 



^g-^W^^-m sm 



in {^z\fuP^^-rrS^ 



dz"(3{z") e*('^+^ 



(46) 



Rewriting the final term by an integration by parts, we get 



:e-*^^"'-"'sm 



in w'^ — m^^ 



1\J Lip- — 

2\/ uj^ — Jo 



dz'p{z')e'''''-'^+'^'^"-"'"^ 
dz'P{z')e"^'^-^'^'"-"''"^ . 



(47) 



Note that the last integral in (|47|) is what we found in ([39]) . where second derivatives were not considered. 

The integration constants can now be determined by the boundary conditions, and in order to do this, let us write 
the complete solution for the axion amplitude. 



(00 + i^O^^'^ 



Substituting (|T7)) and writing explicitly the left and right movers, we have 



(48) 



iz,t) 



_ i{z\/ Lo'^ —ni^ —ujt) 



0O + C + 



c 



aQUO 



2\J Up' — 2\J Up — wi? Jo 

C aou) 



fl\J ijp — m? 2\J up — w? Jo 
Since we do not have incident (or right moving) axions at z = 0, the following condition must be satisfied 

C 



00 + C + 



2\J up — rn? 







while for z > L, there must be only right moving axions, and therefore 



C + aouj / d2'/3(z')e''^''"+^"'""'^ = 



(49) 



(50) 



(51) 



Solving (|5D|) and ([5T|l for the integration constants C, C, and substituting back into (|T7)l . we obtain finally 



(z,t) 



4) aouj 



2^ Up — W? Jo 

aow 



(52) 
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in agreement with the transmitted and reflected waves obtained in ((24|) by the Green function method. 

When the axion mass is very small, so that m << w, and (3{z) is smooth, the reflected axion wave amplitude 
is much smaller than the transmitted axion wave amplitude, because of the rapid oscillation of the phase factor 
gzz (uj+^/uj'^-m'^) ^ This small amplitude is what is neglected in making the WKB approximation of neglecting second 
derivatives of the phases 9, going from to ([55)1 . 



V. ALL ORDERS CALCULATION: DETAILS 



We turn in this section to an all orders calculation of photon-axion and axion-photon conversion, for the case of a 
magnetic field that vanishes for z < and z > L, and is a constant B, so that gB — P is constant, in the interval 
< z < L. We assume a time dependence e^"^* throughout, so that the equations of motion have solutions of the 
form /(z)e^*'^*, where / denotes the fields (p or a; henceforth in this section, the time-dependence factor e^*'^* will 
not be shown explicitly. 

Referring back to ([8]), the photon and axion waves in the region z < have the form 

0(z) = + 0^6-*'=^ , 

a(z) = a/e^"^^ + aije-^"^ . (53) 

Generalizing ^ to allow incoming photon and axion waves from the right (which will be equated to zero later in the 
calculation), the photon and axion waves in the region z > L have the form 

(f,{z) = ^Te'''' + 4>e''''' , 

a{z) = are'"^' + de'"^' . (54) 

In the region < z < L, the propagation eigenmodes are obtained by solving the coupled equations ([5]), which with 
the assumed e~"^' time dependence take the form 

(— — + ■m^)(j){z) = —iujf3a{z) , 

i-LU^-dl)a{z) = iujP(l){z) . (55) 

Assuming an e'^^ dependence of both (/)(z) and a(z), by substituting 0(z) = (f>e'^^^ , a{z) = ae'^^, we find the coupled 
eigenmode equations 

(-0)2 + i^^ + m^)(/) = -iufia, 

{-uj^+K^)a = iujl3<j) , (56) 

which require K to obey the quartic equation 

i-u;^ + K^ + m^){-u;^ + K^)^u;'f3' . (57) 
Solving this, we find that the fields propagate with four possible values for K, either K = ±fc+ or K ~ ±A;_, with 

(fc+)^ = -2-^(i-yTT^) , {k^f^u'-^[i + VTT^) , (58) 

where 

X = 2Puj/m^ . (59) 

For zero magnetic field, we have x = and then fc+ gives the photon mode, and fc_ gives the axion mode. The most 
general solution in the magnetic field region is 

(l>{z) = $+ e'"'=+ + %e''''- + ^+ 6-*^'=+ + e-"''- , (60) 
a{z) = A+ e'''''^ + A- e'''''- + a+ 6-'"'''+ + a- e-'"''- . (61) 



There are restrictions on the integration constants coming from the equations of motion ([56]); they are 

$+=(5A+, ^+^Sa+, Aq^S%, %^S^a . (62) 
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with 



S= ; . 

Thus, the solutions in the magnetic field region take the form 



(63) 



(64) 
(65) 



For a; — > we see that (5 — *■ 0, so that fc+ — ^ cj and fc_ k, and we recover free photons and axions. An important 
property of the above equations is that when the sign of the square root VT+a? is changed from + to — , the 



modes fc+ and fc_ are interchanged, that is, fc+ <-> fc_, and also 5 is replaced by — ^^i^x^ ^ 



= 1/5, that is 



5 <-> 1/5. This property will be used in Appendix B to show that the S'-matrix elements governing all physical effects 
are independent of the branch choice of the square root (a behavior similar to that found in the study of vacuum 
birefringence eS'ects [13, [HI when one does a wave- matching calculation [12,, ilSj] for the case of a rotating magnetic 
field). 

We now impose continuity conditions on fields and spatial derivatives at z = 0, where (3 changes from zero to 
nonzero, and at z = L, where (3 changes back to zero. 



A. Continuity conditions at z = 

At z = 0, the continuity conditions for the fields are 

<Pi + <l)R = 5 A+ + % + 5a+ + , 
ai + an ^ + 6 + + 6 (Pq , 

while the continuity conditions for first spatial derivatives are 

k{(l>i-(t>R) = 5fc+(y4+-a+) + fc_(4>o -(^o) , 
uj{ai - an) = k+ {A'^ - a^) + 6 k-{^Q - ip^) . 



This set of four equations can be recast in matrix form as 



Ml 



ai 



Ma 



\4 ) 



(66) 
(67) 



(68) 
(69) 



(70) 



where the matrices M\ and M2 can be read off from ()66II69|) . and are given in Appendix A. It is straightforward to 
invert the matrix Mi, and defining M12 = Mj"^M2, (j70|) can be rewritten as 




M12 



\4 ) 



(71) 



with the matrix M12 given in Appendix A. 



B. Continuity conditions sX z = L 

At z = the continuity conditions for the fields are 

5 A+ e'^''+ + %e'^''- + 6a+ e-'^''+ + ip^ e~'^^- = <j)T e'^^ 

+ 5 ipQ e" 



A+ e'^^+ + <5 $0 e 



iLk- 



-ikL 



+ e 



iijjL I ~ —iujL 

aT e + a e 



(72) 
(73) 
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while the continuity conditions for first spatial derivatives are 

6 k+{A+ e'^''+ - a+ e-'^''+) + k^{%e'^''- - ip- e-'^^- 



— ALk_ 



) — uj{aT 



(74) 
(75) 



Again, this set of four equations can be recast in matrix form as 







( <Pt\ 


At 


= M4 




\4 ) 




\ « / 



(76) 



where the matrices M3 and A/4 can be read off from (|72ll75p . and are given in Appendix A. And again, it is straight- 
forward to invert the matrix M3, and defining M34 = M^^M^, ([7S)) can be rewritten as 



(77) 







/0t\ 


A + 
^0 


= M34 




V4 J 




[ a J 



with the matrix M34 given in Appendix A. 

C. Construction of the S matrix 

Combining (|7ip and (|77p . we get a matrix relation between the incoming and outgoing wave amplitudes. 




(78) 



From the matrices M12 and Af34 given in Appendix A, explicit calculation gives the S matrix elements tabulated 
in Appendix B. We use them to solve the equation ([TS)) for the following three cases. 



D. Incident photon only 

We assume a photon of unit amplitude incident from the left of 2; = 0, so that a/ = 1, but no axion incident from 
the left, so that 0/ — 0, and no photon or axion incident from the right of z — L, so that a = = 0. 
With this configuration, the set of four equations that must be solved is 

= Sii(j)T + Sisar , 

4>R — S2l4>T + S23aT , 

1 = Szi4>T + 5'33aT , 

aR = S'4i0T + 5'43aT , (79) 
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from which we obtain the transmitted and reflected axion and photon ampUtudes 0t, fpRi o-t, o-r, 



Sis, 



'S'l3'5'31 


— S'll5'33 


'S'2l5'i3 


— 'S'23'5'll 


813,831 


— S'iiiS'33 




^11 


'S'l3'5'31 


— S'll5'33 


>S'4iS'i3 


— S'435'll 


'S'135'31 


— S'll5'33 



(80) 



Referring to (jlOp . these amphtudes obey the unitarity constraint 

l = |aK|' + |aTp + -(|0flP + |0Tn . (81) 



E. Incident axion only 

We assume next an axion of unit ampUtudc incident from the left of z = 0, so that 0/ = 1, but no photon incident 
from the left, so that aj — Q, and again no photon or axion incident from the right of z = L, so that a — = 0. 
With this configuration, the set of four equations that must be solved is 

1 = 8ii4>T + 8i3aT , 

(/"fl — 82l(j)T + 823a,T , 

= S'3i0T + •S'saftT , 
aR = 8ii<j)T + 843aT , (82) 

from which we obtain the transmitted and reflected axion and photon amplitudes (t)T, (f'Rj o-t, O'R-, 

, —833 



813831 — 


5'll<S'33 


'S'23'5'31 ~ 


'S'21'5'33 


813831 — 


'Sll'S'33 


831 


813831 — 


'Sll'S'33 


843831 — 


'S41'S33 


'S'135'31 — 


'Sll'S'33 



'-'43'J31 — '-'41^33 /oo^ 

0-R - -^—^ g, c ■ ^^"^^ 

'J13'-'31 — '-'ll'J33 

Referring again to pU]) . these amplitudes obey the unitarity constraint 



l = |0fl,|' + |0Tp + ^(|aKp + |aTn . (84) 



F. Incident photon and axion 



We assume finally the more general case in which an axion of amplitude </i/ and a photon of amplitude a/ are both 
incident from the left of z = 0, but again no photon or axion are incident from the right of z = L, so that d = cf) = 0. 
In this case, the set of four equations that must be solved is 

0/ = 'Sll0T + 'SlsflT , 
(pR = 82l4>T + 'S23OT , 

a-i = 83i4>T + 'SssttT , 

aR = SiicjiT + 8A3aT , (85) 
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from which we obtain the transmitted and reflected axion and photon ampHtudes (f>T, <t>R^ ^t, cir, 

S13S31 — S'iiS'33 

, _ (<5'23'S'31 — S2lS33)(j)i + {S21S13 — S23Sii)aj 

»J13'-'31 — »-'ll'J33 
5*31'/)/ - 5*110/ 



513531 — 5115"; 



33 



(5435*31 — 54i5'33)(/)/ + (54i5'i3 — 5435'ii)a/ 

O-R = 5 n 5 c ■ ^°^) 

•J13'-'31 ~ »Jll'-'33 

This is of course just a/ times the ampHtudes of (|80p plus (pi times the amphtudes of (|83p . as expected from hnearity. 
The general unitarity constraint of (jlOp now gives, in addition to the constraints on the 5-matrix elements of (j8ip 
and (|84p. the additional constraint 

— fcRe[(5235'3i — 52i533)(52i5i3 — 5235ii)* — 5335^3] 

+ ajRe[(5435'3i — 54i533)(54i5'i3 — 5435ii)* — 5*11531] . (87) 



VI. ALL ORDERS CALCULATION: DISCUSSION 



We now explore various features of the all orders results obtained in the previous section. We begin by showing 
that to leading order in the parameter x — 2ujj3/ni? introduced in (I59p . we recover the result of the Green function 
and WKB/eikonal calculations for the transmitted axion amplitude. When x << 1, expansion of ([55)1 and gives 



from which we find 



5i3 = '7 r [e'^(-^^^)-l] + 0(x^ 



t — j— -^/^ , 1 2 — Tn~2^ 
~? 2 

5ii = 533 = l + 0(x2 



(89) 



Substituting these into (|80p , we get for the transmitted axion amplitude when a photon of unit amplitude is incident 
on the magnetic field region, 



-,2 



4 ^^^237 

which can be rewritten in the usual form given in ([28 



~- (l.e'Hu:-V:j^)-^ , (90) 



<j)T = - ===^6^2^"' ™ ■> sin(-(cj- V^^2-"^2)) • (91) 

We next examine the relationship between the two amplitudes that enter into light through a wall experiments: 
the amplitude (pT for the transmitted axion, when only a photon of unit amplitude is incident [see ()80[) ]. and the 
amplitude ot for the transmitted photon, when only an axion of unit amplitude is incident [see ([83|) ] . We have 

0T(incident photon) _ 5i3 _ _uj_ iL{u>-k) 



At (incident axion) 5*31 k 



g.i.(c^-fc) , (92) 



where we have made use of (BIO). Hence the ratio of magnitudes is loj \JuP- — to all orders in x, for a piecewise 
constant i?, generalizing the result that we found earlier to leading order. 

Finally we turn to the behavior of the axion and photon transmission amplitudes when lo is very close to threshold, 
where we expect to find that the all orders calculation eliminates the unitarity violation that we noted at leading 
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order. Let us parameterize the kinematic variables of the problem in terms of a dimensionless parameter y, which 
takes the value —1 at threshold ui = m, and at the u; value where fc_ vanishes. According to ([55]) . for the product 
k'^kt we have 

klk^_^iuj^-m^/2f-im^/2fil + x^)^ij\ij^-m^-^^) , (93) 

and so fc_ vanishes at to'^ — m? + 0^ . (For values of uj smaller than this, down to threshold dX u) = m, the wave 
number fc_ becomes imaginary.) We thus parameterize lo as follows, 

Lu"^ ^m^ + {l + y)l3^ , (94) 

and assuming that |/3|/m << 1, and that y is of order unity, we get the following approximations for the kinematic 
quantities of interest. 



uj~m + 0{P^), fc=v/l + y|^|, k+~m + 0{l3^), k^ ~ y^\l3\ , 5~~i(i/m . (95) 

The S matrix elements that are needed can now be approximated as follows, working to leading order in |/3|/m and 
in y/l + y, 

^31 -^{l~e-^^n , 
2 m 

5n ^ l--4^[sin(mL)-mL] , 
2m Vl + y 

S33 - I ■ (96) 
Substituting these into (l80|) . we find that |^tP near threshold is given by 

,0 sin^(imL) , , 

(Msin2(imL) + yTT^) +M^{sm(mL)~mLf 

or rewritten in terms of fc = |/3|\/1 + y , 

I^TP- '^^^^i^ . (98) 

(Msin^(i^i) + ^) +Mi(3in(mL)-mL)^ 

Note that once y >> 1, or equivalently k >> |/3p/m, these equations reduce to 



PtI'^ ^ sin^(^mL) , (99) 

cij^ — m-^ 2 

in agreement with the leading order result (|29l ) evaluated at w ~ m. Thus the leading order result is still valid when 
k/m << 1, but requires the correction coming from the all orders calculation when k/m << |/3p/m^. 
Unitarity requires that (fc/w)|0Tp — {k / m)\(j)T\'^ should smaller from unity; from ([98]) we find 



^ , i ,2 ^ ^sin^(^m£) ^ 1 



m / 2M , \ 4 



— |(?;t| 



where we have used the fact that for positive s and t, 

St _l {s + tf -js-t)' ^ l 

(s + t)2 4 (s+f)2 -4 ■ ^ ' 

Using the results of (|92p and (j98p . the squared magnitude of the photon transmission amplitude, for a unit incident 
axion amplitude, is given near threshold by 

,„ -^sin2(imL) 

larP ^ , ^ . (102) 

(^sin2(i^L) + ^) +Mi(sin(mL)-mL)^ 
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Assuming no axion attenuation in the wall, the overall photon transmission probability P{uj) in a light through 
a wall experiment, very near threshold, is the product of |0tP of ([55]) and laxp of (I102p . In the regime where 
<< k/m << 1, P{ui) is well approximated by 

, , /^^sin^fimL) , , 

2m-^(u) — m) 

Defining as the average of P{lo) over an interval extending from w = m to w = m + A, integration of p03p with 

a lower cutoff of kj^ ~ (3"^ /m? , or equivalently, lul — m ^ /3^/(2m^), we find 

P^-/ dcjPicj)^- diu;~m)Piu;)^-^sm\-mL)\ogi——] . (104) 



In a similar fashion, if we define the emerging axion flux from photon to axion conversion in a single B field region, 
for unit incident photon flux, as F{uj) — (fc/m)|(/)Tp, then in the regime where << k/m << 1, F{ll>) is well 

approximated by 

F(c.) ^ sin2(imL) ^ .f' sin^cimi) . (105) 

mk 2 TO3/2^2(a;-m) 2 

Defining F{uj) as the average of F(uj) over an interval extending from 0^ = 171 to lu — m + A, integration of pOSp gives 
the estimate 

1 r^^,_^,_,_, V2p^ ^^.l 



The formulas p04p and (|106p will be of use for estimating the magnitude of photon-axion and axion-photon conversion 
effects in the threshold region. 



VII. MAGNETIC FIELD PENETRATING THE WALL 



Up to this point we have assumed that the magnetic field is present in a vacuum region, where there is no photon 
absorption. Let us now briefly examine the situation where the magnetic field and the photon absorbing "wall" 
overlap. In this case the photon propagation equation must include a complex dielectric constant n — nn + inj with 
nonzero imaginary part nj, and so the eigenmode equations (|56p are modified to read 

(-w^ + + m'^)(f) = -iuoPa , 

{-n^uj"^ +K'^)a = iujPcj) , (107) 

which now require K to obey the quartic equation 

{-uj^ + K^+m^){--n^uj^+K^)^uj^l3^ . (108) 

For P/m << 1, the two types of eignmodes are a "photon-like" mode, where —n^uj^ + is small, so that 

- n^u;^ + = ^5 ~ — ^ !: ^ , (109) 

— LO- + A ^ + TO^ 777/ + (77"^ — 1)W^ 

and an "axion-like" mode, where —uj'^ + + is small, so that 

-n^LO^ + m? + (772 - l)u;2 ^ ' 

Thus, 

■^^photon - njiUj + injuj , 



,/,„2„^2 + ,- "fi"^ ^ii^ till) 

^ +\/ZJ^^7^|m2 + (772 -l)c^2|2 • i^^^J 



15 



We see that for >> m and n/w > 0, the "photon-hke" mode decays as e '^p^otonz^ ^j^jj CTphoton — njuj, while the 
"axion-like" mode decays much more weakly, as e"'^"^'""^, with 

f^axion — f^photon io-,ioo- [^^^ } 

— Ipw^ 

Putting in typical numbers, ui ^ leV, /? ~ 10~^^eV, and (Tphoton ^ lO^cm^^, one has (Taxion ~ lO^^^cm^^, and hence 
there is negligible decay of the "axion-like" wave over a one meter flight path. 



VIII. SUMMARY AND DISCUSSION 



To summarize, we have given three different methods for calculating photon-axion and axion-photon conversion in 
a magnetic field. To get the lowest order transmission amplitude, the WKB method is clearly the simplest, since it 
involves only the integration of a first order differential equation. To also get the reflected wave in leading order, 
the WKB/eikonal and Green function methods are of about equal complexity. To get an all orders answer, a wave 
matching calculation involving 4x4 matrices is necessary. Although this calculation was done only for piecewise 
constant magnetic fields, by expanding the formulas of Appendix B for infinitesimal L as S{/3, dL) = 1 + dLG{(]), and 
using the group property 

S{(3{z < L + dL),L + dL) = S{(3{z < L),L)S{p{L),dL) ^ S{(3{z < L), L)[l + dLG{(3{L))] , (113) 
one gets a differential equation 

dS(p(z < , ^ ^ , , ^ ^, , ^ 

^ \l ' ' =S{piz<L),L)G{piL)) , (114) 

which can be integrated to give the all orders S for an arbitrary longitudinal magnetic field profile P{z). 

All three calculation methods show that when the axion mass is not set equal to zero, there is an enhancement of 
the photon-axion conversion amplitude near threshold a,t uj — m. The lowest order calculations indicate a threshold 
cusp that violates unitarity for small enough fc, and so the all orders calculation is needed to see how, for k very close 
to threshold, unitarity is restored. Deciding whether this threshold enhancement can be used for new experimental 
searches for axions will require further careful analysis. For example, when u) » m the light through walls probability 
obtained from ((29|) is P{uj » m) = (/3L/2)'*, which for large values of the effective path length mL is much larger 
than the probability (|104p for the same geometry when w is near threshold, assuming that the effective bandwidth A 
is of order the axion mass m. However, further investigation will be needed to see if (|104p suggests ways of exploiting 
the threshold enhancement by using very small laser bandwidths A << m. Similarly, further analysis will be needed 
to determine whether the enhanced axion flux near threshold given by (|106p has astrophysical implications in strong 
magnetic fleld environments. 



IX. ACKNOWLEDGMENTS 



The authors wish to thank E. I. Guendelman for his participation in initial phases of this project, and for helpful 
email correspondence suggesting a WKB approach and pointing out the role of the k factor in the unitarity relation. 
The work of SLA was supported in part by the U. S. Department of Energy under Grant No. DE-FG02-90ER40542, 
JG and FM were partially supported by grants from FONDECYT 1050114 and 1060079. JL-S acknowledges the 
support from MEC/FULBRIGHT-FU2006-0469. 



X. ADDED NOTE 



Carlo Rizzo has raised the pertinent question of how the factor sin^(imL), and its square, in the formulas 
(103)-(106), are to be interpreted when mL » 1. This factor is the evaluation at threshold of the factor 
sin^ (^-^('^ ~ Vw^ — m^)) ~ sin'^ (iL(m — ^2m{u) — m))) . When ^L^/2m/S. » 1, the sine function has many 
oscillations over the integration interval in (106), (104), and it then can be replaced by the respective averages 
(sin^ (ii(m - yj2m{uj-m))))Av = 1/2, and (sin* {\L{m - ^2m{uj - 'm))))Av = 3/8, in (106) and (104) respec- 
tively. When iL\/2mA is of order unity, the full w-dependent form of the argument of the sine function should be 
included in the integrals in (106), (104) in place of its threshold evaluation, and the integrals can then be computed 
numerically. 
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APPENDIX A: MATRICES FOR MATCHING AT ^; = AND z 



The matrices Mi, M2, and M12 = ^M^ needed for the 2; = match are 



Ml 



Mo 



Ml 



12 



110 
11 
fc -fc 
w -w 



1 

V (5 fc_ 



1 

-5 k 



1 



1 

A;_/fc 1 



<5 
1 

-(5 fc+ 

- k^jk 



1 [ l-k-/k l + k-/k 

2 5(l + A;_/a;) (5(l-fc_/a;) 
\^(l-fc_/a;) 5ll + k-/uj) 



5{l + k+/k) 5{l-k+/k) 

5{l-k+/k) 5{l + k+/k) 

l + k+/io l — k+fu) 

l — k+/u l + k+/uj 



(Al) 



The matrices M3, M4, and M34 



M3 ^M4 needed for the z 



L match are 



Ms 



M4 



M34 





/ 6*^'=- e" 


-iLk— 








— iLk- 






fc_ e'^'=- -k- 


^—iLk- 






\5 k- e'^''- -5 k. 


^—iLk- 






( ^iLk ^—iLk 












^iLu) 
















y 







jg-iLfe+ 
g-iLfc+ 

-k+ 6-'^*=+ 



1 



2(1 -(52) 
/ {l + k/k. 



iL(k-k-) 

(l-fc/fc_)e'-f'(*^+'^--) 
-5(l + A;/fc+)e'^^('=^-'=+) 
-5(l-fc/fc+)e'^('=+'=+) 



(3 



-iL(fc+fc_) 

iL(k--k) 



(1 - 

(l + fc/A:_ 
-5(l-A;/fc+)e-'^^('=+'=+) 
-5(l + A;/A;+)e'^('=+-'=) 



-Sil+Lo/k-)e'^^-^-'''-^ 
-8{l-uj/kJ)e'^^'^+^-^ 
(1 +w/fc+)e*-^('^-'=^+) 
(1 - a;/fc+)e'^('^+'=+) 



-(5(l-w/fc_)e-*^('^+'=-) \ 
-5(1 +w/fc_)e'-^(*=-~'^) 
(1 - w/A:+)e-*-^("+'=+) 
(l+a;/fc+)e'^('=+-'^) / 



(A2) 
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APPENDIX B: MATRIX ELEMENTS OF 5i, 



2 {5^ - 1) k_k+k Sn = e*" 



2 {5^ - 1) k-k+k Si2 
2(6^^1) k-k+k Si3 



521 



2{S^ - 1) k-k+k S2 



2(6"^ -l)k- k+kS; 



23 



•Ssi 

<S'32 



2 - 1) cok-k+S: 



2{6^ - l)wfc-A:+534 



'41 



'S'42 

5'43 

2 - 1) iok-k+Si4. 



le 



ik+{kt +k'^)sm{Lk-) - 2k+k-kcos{Lk-) 
(5^ [ik-{k+ + fc^) sin(ifc+) - 2A;+ A;_ fccos(Lfc+; 

- sin(ifc_) + 5^ k-ik"^ - fc^) sin(Lfc+) 
(w + k)k-k+[cos{Lk-) — cos(LA:+)] 
i [k+ {k^ + Lok) sm{Lk- ) - k- (fc^ + tok) sm{Lk+)] 



-iLk 



2 {d^ - 1) k-k+k Si4 = 6e 



— iLu) 



{uj — k)k-k^[cos{Lk^) — cos{Lk- 



—i[k-\-{k'^ — Luk) sm{Lk-) — k-{k'^ — u)k) sin(LA;+) 



-e bi2 



-iLk 



k+[2k-kcos{Lk-) + i{k'^ +fc^)sin(Lfc_ 
+S'^k-[2k+kcos{Lk+) + i{kl + k^) siii(ifc+)] 



k+[k- {u! — k) cos (Lk-) — i {h?L — cuk) sin (Lk-)] 



-k- [k+ (w — k) cos {Lk+) — i [k+ — Lok) sin {Lk+)] 



2 {S^ - 1) k-k+k S24 = e 



-iLui i 



k-[k+{uj + k) cos{Lk+) + i{k+ +Lvk) sin(iyfc+)] 



+k+[k-{w + k) cos{Lk-) + i{k'i + wA;)sin(Lfc_)] 



-iL{uj-~k) 



k 



'13 



33 = «e 



g-il,(a;+fe) ^ g^^ 

k-[2iujk+ cos{Lk+) + (w^ + k+) sin(Lfc+)] 
-5"^ k+ [2iu> k- cos(Lfc_ ) + (w^ + fc^ ) sin(Lfc_ )] 



-ILlj 



k-{k\ - a;2) sin(Lfc+) - 5^ k+{kl. - w^) sin(LA:_) 
k 



_^iH.-k) k 

LO 

—e D34 



-iLuj 



k- [2a; k+ cos(Lk+) + i{u>'^ + k\) sin(LA;+) 
+5'^k+[2oj k- cos{Lk-) + i{iJ^ + kV) sin(Lfc_)] 



(Bl) 
(B2) 

(B3) 

(B4) 
(B5) 

(B6) 
(B7) 

(B8) 

(B9) 
(BIO) 

(Bll) 

(B12) 

(B13) 

(B14) 
(B15) 

(B16) 
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The matrix elements of S have a large number of symmetries that arise from the structure of M12 and M34. First, all of 
the Sij are even functions of A:+ and even functions of fc_. Second, the matrix elements above have one of the following 
two structures: (i) [S/{l~S'^)][f{k+,k^,k,Lu)^f{k^,k+,k,uj)], (ii) - S^)][S^ f{k+, k_, k, to) ~ fik^, k+, k, to)], 

with the functions / different for each matrix element. Since we have seen that under the interchange A;+ <-> /c_ one 
has S <-> 1/(5, these two structures are both invariant under the interchange of fc-|- and fc_, and so all matrix elements 
Sij are invariant under this interchange. Finally, there are a number of relations between matrix elements under 
reversal of the sign of fc or w, as follows 

Sll = S22\k^-k , S21 — Si2\k^-k , S32 = <5'3l|fc^-fe , 5*13 = 5'23|fc^-fe , 5*41 = 5'42|fe^-fe , <S'l4 = <5'24|fc^-fe , 
S33 = 844]^^-^ , S'34 = 843]^^-^ , S'13 = 814]^^-^ , 523 — S24\uj^-uj , S41 — S'siloJ^-w ! S42 = 5*321^^-0) 

(B17) 

These relations serve as useful checks on the calculation. 
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FIG. 1: Light shining through a wall setup. The wavy line indicates the photon field and the dashed line the axion field. In 
region I, there is only the photon wave, in region II the magnetic field gives rise to an axion wave, which exits with the photon 
wave into region III. In region IV, the "wall", the photon wave is absorbed, leaving only the axion wave in region V. In region 
VI, the magnetic field regenerates a photon wave from the axion, which exits into region VII, where the regenerated photon is 
detected. 



